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Abstract. We implement a Monte Carlo radiative transfer method, that uses a large number of monochromatic luminosity 
packets to represent the radiation transported through the system. These packets are injected into the system and interact 
stochastically with it. We test our code against various benchmark calculations and determine the number of packets required 
to obtain accurate results under different circumstances. We then use this method to study cores that are directly exposed to the 
interstellar radiation field (non-embedded cores). Our code predicts temperature and intensity profiles inside these cores which 
are in good agreement with previous studies using different radiative transfer methods. 

We also explore a large number of models of cores that are embedded in the centre of a molecular cloud. We study cores with 
different density profiles embedded in molecular clouds with various optical extinctions and we calculate temperature profiles, 
SEDs and intensity profiles. Our study indicates that the temperature profiles in embedded cores are less steep than those in 
non-embedded cores. Deeply embedded cores (ambient cloud with visual extinction larger than 15-25) are almost isothermal at 
around 7-8 K. The temperature inside cores surrounded by an ambient cloud of even moderate thickness (Ay ~ 5) is less than 
12 K, which is lower than previous studies have assumed. Thus, previous mass calculations of embedded cores (for example in 
the p Ophiuchi protocluster), based on mm continuum observations, may underestimate core masses by up to a factor of 2. 
Our study shows that the best wavelength region to observe embedded cores is between 400 and 500 yum, where the core is 
quite distinct from the background. We also predict that very sensitive observations (~ 1 - 3 MJy sr"') at 170-200 fjm can be 
used to estimate how deeply a core is embedded in its parent molecular cloud. The upcoming Herschel mission (ESA, 2007) 
will, in principle, be able to detect these features and test our models. 
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1. Introduction 

The initial stages of star formation are not very well under- 
stood. The general view is that molecular cloud cores with sizes 
~ 0.1 pc and masses of a few Mq, collapse to form stars, either 
in isolation or in clusters. Various observations have been asso- 
ciated with different stages of this scenario (see Andre, Ward- 
Thompson & Barsony 2000). Class objects coiTespond to 
the first stage in the evolution of a protostar, where a central 
luminosity source has been formed in the centre of the core, 
and matter accretes on to it. Class I objects correspond to a 
later stage of collapse where a disc has started to form around 
the central object but there is also a residual surrounding en- 
velope. Accretion onto the central protostar continues but at 
a lower rate. Class II and Class III objects correspond to the 
classical T Tauri (CTT) and weak-line T Tauri (WTT) stars, re- 
spectively. CTT stars have well defined discs, whereas in WTT 
stars the inner discs have dissipated. In addition, observations 
from Earth (IRAM, SCUBA/JCMT) and space {IRAS, ISO) of 
various molecular clouds have revealed condensations that ap- 
pear to be on the verge of collapse or already collapsing (e.g. 
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Myers and Benson 1983, Ward-Thompson et al. 1994, Ward- 
Thompson et al. 2002, Kirk 2002). These condensations are 
referred to as prestellar cores. 

Isolated prestellar cores have been observed inside molec- 
ular clouds (e.g. L1544; Ward-Thompson et al. 1999). These 
cores are considered to be precursors of isolated low mass star 
formation. Isolated prestellar cores have extent t 1.5 x 10^ AU 
and masses 0.5 - 35 Mq (Andre, Ward-Thompson & Barsony 
2000). They are not in general spherically symmetric and they 
appear to have flat central density profiles. Magnetic fields are 
also present and they may play a role in core stability (Andre, 
Ward-Thompson & Barsony 2000 and references therein). 

Prestellar cores have also been observed in young proto- 
clusters, such as p Ophiuchi (Motte et al. 1998, Johnstone et al. 

2000) and NGC 2068/2071 (Motte et al. 2001, Johnstone et al. 

2001) . p Oph is a star-forming cluster of about 1 pc diameter, 
with estimated average particle density n(H) ~ 2-4x 10"* cm"^ 
and thermal gas pressure ~ 10* cm"^ K (Liseau et al. 1999). 
In this region there have been detected 100 structures, 59 of 
which are identified as prestellar cores and the remaining as 
embedded young stellar objects (Motte et al. 1998). The extent 
of the prestellar condensations is 2 - 4 x 1 0^ AU (more compact 
than isolated prestellar cores), and they have sharp edges. Their 
estimated masses are 0.05 - 3 Mq. NGC 2068/2071 are proto- 
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clusters in the Orion B cloud complex. Observations (Motte 
et al. 2001, Johnstone et al. 2001) have revealed a filamentary 
structure with ~ 70 starless condensations having small sizes 
(~ 5000 AU) and masses from ~ 0.3 Mq to ~ 5 Mq. 

Mass estimates from mm continuum observations, where 
the cores are optically thin, suggest that the initial mass func- 
tion (IMF) could be determined by fragmentation at the pre- 
stellar stage of star formation (e.g. Andre, Ward-Thompson 
& Barsony 2000). The question of whether fragmentation 
can produce the smallest masses in the IMF is still open. 
Observations of very low mass prestellar condensations are 
crucial for answering this question but they are beyond the lim- 
its of today's telescopes. Furthermore, current mass estimates 
are uncertain, due to our limited knowledge of the properties of 
the dust in and around these cores, and of the dust temperature. 
Previous studies have assumed isothermal dust at 12-20 K (e.g. 
Motte et al. 1998, Johnstone et al. 2000). More recent radiative 
transfer studies (Evans et al. 2001, Zucconi et al. 2001) model 
cores that are illuminated directly by the isotropic interstellar 
radiation field and find that the temperature decreases towards 
the centre of the core. However, these studies cannot be applied 
to embedded cores, because in this case the illuminating radi- 
ation field is not the interstellar one, and, in general, it is not 
isotropic. 

In this paper, we present a more realistic model that treats 
cores that are embedded in molecular clouds. We use a Monte 
Carlo radiation code we have developed to study cores approx- 
imated by Bonnor-Ebert (BE) spheres. In Section|2 we discuss 
the basics of our code and the tests we have performed to check 
its validity. In Section |3 we discuss how we adapt our code 
to treat the radiation transfer in externally illuminated spheres, 
and present the tests we have performed. In Section B31 we 
briefly examine the effect of different dust properties on our 
calculations and, in Section |3 we study BE spheres exposed 
directly to the Black ( 1 994) interstellar radiation field and com- 
pare our results, which were acquired by a different radiative 
transfer method, with those of Evans et al. (2001) and Zucconi 
et al. (2001). In Section |5] we study the more realistic case 
of cores embedded in molecular clouds; we calculate the dust 
temperature in their interiors, their spectra and their intensity 
profiles at different observing wavelengths. Finally, we sum- 
marise our results in Section|6] 

2. Monte Carlo Radiative Transfer 

We implement a method for radiation transfer calculations 
based on a Monte Carlo approach, similar to that developed 
by Wolf et al. (1999) and Bjorkman & Wood (2001). We make 
use of the fundamental principle of Monte Carlo methods, ac- 
cording to which we can sample a quantity ^ e [^1,^2]. from a 
probability distribution using uniformly distributed random 
numbers K e [0, 1], by picking ^ such that 



(1) 
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Here we briefly outline the basics of our code, named 
Phaethon, after a Greek mythical hero. 



We discretise the radiation field of a luminosity source (star 
or background radiation) using a large number of monochro- 
matic luminosity packets (hereafter referred to as 'L-packets'). 
The frequency of an L-packet is chosen from the source ra- 
diation field ly, using Eq. [2 which, in this case, becomes 

Iy(T)dv/£' Iy(T)dv = "Ry, (K ^ [0,1]). Each of the L- 
packets is injected stochastically into the medium, either from 
a specific point (for a point star) or from the boundaries of 
the system (for background radiation). For an isotropic radia- 
tion field from a point source (e.g. Yusef-Zadeh et al. 1984), 
the direction (0,<p) of the L-packet is calculated using 6 - 
cos-i(l - 2??fl) and = 27fR^ {Kg ,1^^ e [0, 1]). Each L- 
packet is also assigned an optical depth, using Ty - - In Hr, 
Kj e [0, 1], and this determines how far the packet propagates 
into the the medium before it interacts with it. 

The computational domain in which the L-packets propa- 
gate is divided into a number of cells. In regions where the 
density or the temperature gradients are large, more cells are 
needed. We can fulfil both conditions by constructing cells with 
dimensions 5 cell less than, or on the order of, the local direc- 
tional scale-heights. 



S, , II i Mm [hp, hr} 
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In the direction given by the unit vector k, the directional scale 
heights are 
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In theory, we can construct a grid with a very large number of 
cells to satisfy our accuracy requirements in calculating tem- 
perature. However, if we use a large number of cells we need 
a large number of L-packets to interact with each of these cells 
so that the statistical noise of our calculations (of the order of 
1 / y/Nabs, where A^abs is the number of L-packets absorbed in 
each cell) is small, that will increase the computational time of 
our calculations. 

If Ttotai is an L-packet's total optical depth then in order to 
calculate the distance it propagates into the system before it 
interacts with it, we need to calculate the line integral along the 
path of the packet. 
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In the general case it is not possible to calculate the preceding 
integral analytically. Our approach is to approximate this inte- 
gral with a sum: AS - ^APi) - Tji^Si. The element 
step 6S i that each L-packet propagates should be small, so that 
the density remains almost constant along this step. Also, the 
element optical depth, 5t, - k,\ pi SS ,, should not be larger than 
the remaining total optical depth of the L-packet r,; t, is just 
the optical depth that the L-packet still has to propagate after / 
steps, Tj = Ttotai - 2)=o ^''"i ■ To satisfy the above requirements 
we chose an element step according to the following condition: 



SS i = MIN [j]p hp, T] I, (t,- + e) 1, 77, 



(5) 



where / = {kaPi) and rjp, 77, rjr are constants that determine 
the accuracy we demand (typical values are between 0. 1 and 1). 



D. Stamatellos and A. P. Whitworth: Monte Carlo Radiative Transfer in Embedded Prestellar Cores 



3 



Initialize grid 



Calculate 
Dust Opacities 



Launch L-packet 

Calculate L-packet frequency 
r Calculate L-packet direction ^ 

( Calculate L-pacto optical deplli ^ 



Propagate L-packet 



If L-packet insidey If L-packet outside 



the system 



the system 



albedo<R albedo>R 



Absorb and re-emit 
L-packet 

Calculate cell new temperature^ 
^ Calculate L-packet frequency j 
f Calculate L-packet direction ] 



fcalculflte 



L-packet optical depth 



D 



Scatter L-packet 
f Calculate L-packet direction 

\ Calculate L-packet optical depth j 

Re-remit L-packet 



T 



Place into frequency 
and direction bins 



Fig. 1. phaethon: Code flow chart. Luminosity packets are in- 
jected into the system, propagate, interact and finally escape. 
The propagation routine is the most computationally expensive 
routine. 



The first term (jjp hp) ensures that the density does not change 
much in one element step {hp is the density scale height in the 
direction that the L-packet propagates), the second term (// 1) 
ensures that the element step is less than the mean free path of 
the L-packet and the third term [(t, + e) /] takes effect on the last 
step (e is a very small number). The last term (rfr |r|) ensures 
that the distance the L-packet travels in one element step is less 
than the distance from the luminosity source. This term comes 
into effect when a gap exists around the source. The smaller the 
factors Tjp, T], rjr are chosen, the better the accuracy in propagat- 
ing L-packets , but on the other hand a smaller element step 
means more computation. We propagate the L-packet follow- 
ing the above procedure until t,- < or until the packet escapes 
from the system. 

When an L-packet reaches an interaction point within the 
medium (at the end of its optical depth), it is either scattered or 
absorbed, depending on the albedo. If the packet is absorbed its 
energy is added to the medium and raises the local temperature. 
To ensure radiative equilibrium the L-packet is immediately re- 
emitted. The new temperature of the cell that absorbs the packet 
is found by equating the absorbed and emitted energies. The 
re-emitted L-packet has the same energy but a new frequency 
chosen from the difference of the local medium emissivity be- 
fore and after the absorption of the packet (Bjorkman & Wood 
2001). The direction of the reemitted L-packet is random. If 
the L-packet is scattered then it is assigned a new direction, us- 
ing the scattering phase function due to Henyey & Greenstein 
(1941). Then the packet propagates again in the medium to a 
new interaction location. This procedure continues until all the 
packets escape from the system. They are then placed into fre- 
quency and direction-of-observationbins. 

The general flowchart of the radiative transfer code 
Phaethon is shown in Fig.^ The shaded parts in this diagram 
refer to procedures that, in general, need to be done more than 
once for each L-packet and consequently are those that dictate 



the efficiency of the code. A simple code efficiency analysis, 
indicates that the L-packet propagation routine takes about 25- 
50% of the computational time, depending on the specific prob- 
lem. Thus, this is the routine that should be targeted by any ef- 
forts to diminish the running time of the code. Time efficiency 
is very important since a large number of L-packets is needed 
for good statistical results. To reduce the running time of the 
code, whilst maintaining good results, the specific nature of the 
system we examine should be taken into account (for exam- 
ple, in the case of a uniform density sphere we can propagate 
each L-packet in a single step), any kind of symmetry should be 
exploited (e.g for spherically symmetric systems) and look-up 
tables should be used to solve for the cell temperature after ab- 
sorbing a packet and then to calculate the reemission frequency 
of the new packet. 

We have tested our code against benchmark calculations 
proposed by Ivezic et al. (1997) for a star surrounded by a 
spherical envelope, first with constant density, and then for 
density decreasing as r"^. They used three different, well es- 
tablished radiative transfer codes using different numerical 
schemes to solve a set of benchmark spherical geometry prob- 
lems. In all cases, these methods gave differences smaller that 
0.1% and, as Ivezic et al. noted, the solution should be consid- 
ered exact. Our code reproduces those results and also the re- 
sults of Bjorkman & Wood (2001), for a disk-like structure em- 
bedded in an envelope. These tests demonstrate the validity of 
our radiative transfer code. We will not present these tests but, 
instead we will later discuss two different tests: (i) the 'thermo- 
dynamic equilibrium test' in which a Bonnor-Ebert sphere is 
illuminated by a uniform blackbody radiation field and (ii) the 
'pure scattering test' in which the albedo of the dust is set to 1 . 

This Monte Carlo radiative transfer method conserves en- 
ergy exactly, accounts for the diffuse radiation field, can be 
implemented for any geometrical structure and is very effi- 
cient, making it very attractive for use in a variety of problems. 
However, it can be implemented, without iteration, only when 
the opacity is independent of temperature, so the method is use- 
ful for treating radiation transport against opacity due to dust 
grains which are large enough to be in thermal equilibrium. 

3. Radiative Transfer in Prestellar Cores 

3.1. Core Density Profile 

A simple approach to prestellar cores is to assume that they are 
isothermal spheres in which gravity is balanced by gas pressure 
(Bonnor-Ebert spheres; Bonnor 1956, Ebert 1955). Recent ob- 
servations (e.g. Alves et al. 2001, Ward-Thompson et al. 2002) 
show that this is a good approximation for many cores. We use 
the Monte Carlo radiative transfer code to study cores embed- 
ded in an isotropic interstellar radiation field. We choose to 
parameterise BE spheres using the temperature, the mass of 
the sphere, and the external ambient pressure on the sphere. 
This type of parameterisation is quite useful when examining 
prestellar cores in the same molecular cloud, in as much as we 
can presume that they all experience roughly the same exter- 
nal pressure. The sphere is divided into a number of concentric 
cells (typically 50-100) with equal radial width (Fig. IS. 
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Fig. 2. Schematic view of Bonnor-Ebert sphere model.L- 
packets are ejected from the point (0, 0, R) at such an angle as 
to imitate an isotropic radiation field. 

In our study we assume isothermal gas BE spheres. 
Generally, dust and gas do not have the same temperature un- 
less the density is quite high, in which case they are thermally 
coupled (n > 1 - 3 X 10"* cm"^ Mathis et al. 1983; Whitworth, 
Francis and Boffin 1998). However, even non-isothermal mod- 
els that allow a small gas temperature gradient, give density 
profiles that are very close to the BE profile (Evans et al. 2001). 
Thus, our results for isothermal spheres should represent non- 
isothermal spheres reasonably well. 

3.2. The Illuminating Radiation Field 

Because the core is spherically symmetric and the radiation 
field is isotropic, we can -without loss of generality- inject all 
L-packets at the point (0, 0, R) along the yz plane, where R is 
the radius of the BE sphere (Fig. If A) is the integrated in- 
tensity of the radiation field, then the total luminosity incident 
on the sphere is Ltotai = ^lo 47r/?^ . If we use A^^ luminosity 
packets, the luminosity per L-packet is 6L = (jtIq AnR^) j N p . 
For isotropic intensity the injection angle probability is 

ped0 = 2 cos {&) sin (e)de , ^ < < n , (6) 

and the L-packet injection angle is therefore 

e = cos-' [-'Rg"^], He e [0, 1]. (7) 

For the spectrum of the radiation incident on the core, 
we use the Black (1994) interstellar radiation field (hereafter 
BISRF). Black has compiled an average Galactic background 
spectrum from radio frequencies to the Lyman continuum limit, 
based both on observations and theoretical modelling (Fig.|3j. 
This spectrum consists of an optical component with a peak 
at around 1 /vm, due to radiation from giant stars and dwarfs; 
a component due to thermal emission from dust grains with 
a peak at around 100 jum; mid-infrared radiation from non- 
thermally heated grains in the range 5-100 jum; and the cos- 
mic background radiation with a peak around 1mm {T - 
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Fig. 3. Black (1994) Interstellar Radiation Field. 

2.728 ± 0.004 K). This background is similar to that of Mathis 
et al. (1983) apart from the region from 5 to 400 yum, where 
it is stronger on the basis of COBE data. As noted by Black, 
his estimate only accounts for continuum radiation and does 
not include strong emission lines, which may have significant 
power in the far-IR and submillimetre part of the spectrum. 

The BISRF seems to be a good approximation to the inter- 
stellar radiation field in the solar neighbourhood. However, it 
is not always an appropriate choice when studying prestellar 
cores, because many cores are embedded in molecular clouds. 
Consequently, the radiation field is attenuated at short wave- 
lengths (< 30 - 40 fim) because the surrounding cloud ab- 
sorbs a large part of this radiation, and enhanced at long wave- 
lengths (> 50 jum) due to the thermal emission from the molec- 
ular cloud (Mathis et al. 1983). Also the radiation field may be 
anisotropic. In this work, initially we study cores directly ex- 
posed to the BISRF (like the previous studies of Evans et al. 
2001 and Zucconi et al. 2001) but we also extend our study 
to the more realistic case of cores inside molecular clouds of 
diff'erent sizes. 

3.3. Dust Opacities 

Typical dust temperatures in prestellar cores are quite low (5 - 
20 K) and under these conditions dust grains are expected to 
coagulate and accrete ice mantles. Following recent studies of 
prestellar cores (Evans et al. 2001, Zucconi et al. 2001), we 
use absorption opacities calculated by Ossenkopf and Henning 
(1994) (hereafter OH) for a standard MRN (Mathis, Rumpl & 
Nordsieck 1977) interstellar grain mixture (53% silicate & 47% 
graphite), with grains that have coagulated and accreted thin 
ice mantles over a period of 10^ years at densities 10''cm-''. We 
also assume a gas-to-dust mass ratio of 100. 

Ossenkopf and Henning only calculated absorption opaci- 
ties down to 1 yum, so below this value we use the MRN stan- 
dard model (grains without ice mantles) with optical constants 
from Draine & Lee (1984), after scaling to fit the OH absorp- 
tion opacity at 1 yum (Fig.0}. In any case, the choice of absorp- 
tion opacities below 1 fim, does not play an important role in 
our calculations since at these short wavelengths the core is op- 
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A(/^m) 

Fig. 4. Ossenkopf & Henning (1994) + MRN (1977) opacities. 
The solid line represents the absorption opacity and the dashed 
line the scattering opacity (/Cabs + A^scat = 10.45 x 10^ cm^g"' 
at /( = 0.55 yum). 



tically thick and the radiation does not penetrate much inside 
the core. 

Also due to lack of data for scattering opacities we use 
the MRN scattering opacities after scaling them as before. We 
should note though that for dust in prestellar cores, scattering 
is expected to be less by at least a factor of 2 (Ossenkopf; pri- 
vate communication). The choice of scattering opacities does 
not greatly affect the temperature in the inner regions of the 
core, since scattering is only important for short wavelength 
(520 yum) radiation, which, anyway, cannot penetrate deep in- 
side the core. We will discuss the effect of scattering in more 
detail later. 
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Fig. 5. Thermodynamic equilibrium Test for a Bonnor-Ebert 
sphere, (a) Temperature versus distance from the centre of the 
sphere, (b) Spectrum of the incident and the emergent radiation 
in units of B = ctT'^/tt (there is no difference at all). 
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Fig. 6. (a) SED of the incident and the emergent radiation (no 
difference at all), (b) Pure scattering test: Intensity profiles at 
wavelengths 850, 200, 450 and 0.55 /im (top to bottom), for 
mean scattering cosine g = 0.5 . The dotted lines against the 
righthand margin of the right plot correspond to the background 
radiation at each of the above wavelengths. For A - 450 i-im the 
statistical noise is larger because fewer L-packets are emitted at 
this wavelength than at other wavelengths. 



3.4. Code Tests 

3.4.1 . Test 1 : Thermodynamic Equilibrium 

Consider a system that is illuminated by a uniform, isotropic 
blackbody radiation field of temperature T. Thermodynamic 
equilibrium dictates that every part of the system will adopt 
the same temperature T. This also means that the intensity of 
the radiation coming from the system is the same as that of the 
illuminating blackbody field. It is easy to see this from a simple 
radiative transfer calculation. If /v(0) = By{T) is the intensity 
of the incident radiation at a specific direction on the system, 
the intensity Iv(D) of the radiation that escapes, after travelling 
distance D inside the system, is 

IAD) = /v(0)e-^'(^> + BAT)[l - e-^'<^>] = By . (8) 

Practically this means that the system is invisible to an ob- 
server This test can be applied to any structure (e.g. spheres, 
discs, non-symmetric structures) and it is a good way to check 
the main radiative transfer code (i.e. L-packet injection, prop- 
agation, absorption, temperature correction and reemission). It 
is a very discriminating test and we suggest that it should be 
apphed to all radiative transfer codes. 



We perform the thermodynamic equilibrium test for an un- 
stable Bonnor-Ebert sphere (^out=11.8, M - 4.5 Mq, T - 
11 K, Pext = lO'* cm"^ K). Initially, we do the test with a 
blackbody illuminating field having T = 10 K and then with 
r = 20 K (using lO' luminosity packets). As seen in Fig. |5] 
the output spectrum is the same as that of the illuminating field 
and the temperature at any distance from the centre of the core 
is constant and equal to that of the radiation field. Small vari- 
ations on the order of 0. 1 K are not important and are due to 
statistical noise. 

3.4.2. Test 2: Pure Scattering 

If radiation field incident on the sphere is isotropic and if the 
L-packets just pass through the sphere without interacting, the 
observed intensity will be the same at each impact parameter 
b, and equal to the intensity of the illuminating field. The same 
holds if the L-packets just get scattered, i.e. when the albedo of 
the grains is set equal to 1 . It is easy to understand this when 
the scattering is isotropic, but the same is also true for non- 
isotropic scattering. Since the incident field is isotropic, the L- 
packets come from aU directions and the eff'ect of scattering 
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Table 1. Non-Embedded Prestellar Cores: Model Parameters 



model ID 




T''( a:) 


M-^CMq) 


Sout 


' V 


BEl 




10 


2 


2.6 


4.0 


BE2 




10 


3.5 


4.5 


8.3 


BE2.2 


l(f 


10 


3.5 


9.9 


24.0 


BE3 


10'^ 


15 


2 


1.7 


2.4 


BE4 


l(f 


15 


4 


2.4 


3.6 


BE5 


l(f 


15 


6 


3.2 


5.2 


BE5.2 


10* 


15 


6 


21.6 


64.0 


BE6 


10' 


15 


1 


2.1 


9.8 


BE7 


10' 


15 


2.6 


5.0 


30.1 


BE7.2 


10' 


15 


2.6 


8.5 


62.0 


BE8 


5 X 10* 


15 


2 


2.6 


8.9 


BE9 


5 X 10'' 


15 


3.5 


4.5 


18.2 


BE9.2 


5x lO'' 


15 


3.5 


10.1 


56.2 



" External pressure 
Gas temperature 
Bonnor-Ebert sphere mass 
^ parameter (sphere is supercritical if f > 6.451) 
Visual optical depth to the centre of the sphere. 



will simply be to rotate the whole radiation field through an 
angle but the field will remain isotropic. The same argument 
holds if an L-packet undergoes more than one scattering. Thus, 
if the radiation just gets scattered in the medium, the emergent 
spectrum will again be the same as that of the incident radia- 
tion. 

We perform this test for a BE sphere (parameters: ^out = 
4.1, M = 4 M0, r = 11 K, Pext = 10^ cm-^^ K). This time 
the sphere is illuminated by the BISRF. We do calculations for 
mean scattering cosine (isotropic scattering), 0.5 and 1 (us- 
ing 2 X 10^ luminosity packets). We present our results for the 
g = 0.5 case in Fig. |6l The code successfully passed this test 
too. 

4. Non-Embedded Prestellar Cores 

We use the term non-embedded piesteWai cores to refer to cores 
that are directly exposed to the BISRF. We perform simulations 
for a number of Bonnor-Ebert spheres under different exter- 
nal pressures and for various gas temperatures and masses. In 
Tabled we list the parameters of our models to show the pa- 
rameter space investigated. BE spheres with the same set of 
Pen, T and M correspond to one subcritical and one super- 
critical sphere. These can be distinguished by the ^out value; 
if ^out > 6.451 then the sphere is supercritical. For each model 
we calculate the temperature profile of the dust in the core, the 
core SED and intensity profiles at different wavelengths (90, 
170 and 450 //m). 

4.1. Temperature Profiles 

The dust temperature inside the core drops from around 17 K 
at the edge to a minimum at the centre, which may be as low 
as 7 K, depending on the total optical depth of the sphere. The 
higher the optical depth to the centre of the core the lower the 
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Fig. 9. Components of the spectrum for the BE2.2 model. The 
dotted black line is the SED of the radiation incident on the 
core, the solid line is the output SED, the dash-dot line is the 
part of the radiation that passes through the core without inter- 
acting, the short-dashed line is the core emission and the long- 
dashed line is the scattered light. 

central temperature and the larger the temperature gradient. In 
Fig. we plot dust temperature profiles for three representa- 
tive core models with different density profiles. 

Our results are in general agreement with previous similar 
calculations by Evans et al. (2001) and Zucconi et al. (2001). 
We compared our results with those of Evans et al. for a sys- 
tem as close we could to get to one of their models. They do 
not mention what opacities they use for A < I /im, they use an 
ISRF at A < 1, different from the BISRF, and some parame- 
ters in their model are unclear We find that the temperature we 
calculate at the edge of the core is 3 K higher than their calcu- 
lation (17 K rather than 14 K). This difference can be explained 
in terms of different opacities and different ISRF for A < I fj.m. 
We also find that the temperature is almost 1 K lower than re- 
ported by Evans et al. at the centre of the core. This may in part 
be due to slightly different density profiles. Zucconi et al. also 
reported a higher estimated temperature at the edge (» 16.5 K) 
and a lower temperature (^ 0.3 K lower) at the centre of the 
core, when they compared their model with the Evans et al. 
calculations. However, these are small differences. 

4.2. SEDs and Intensity Profiles 

We see from the component version of the SED (Fig.|9}, where 
we plot the contribution from scattered, processed and direct 
photons to the SED for the BE2.2 model, that short wavelength 
radiation (A ; 50 jum) is absorbed from the core and then is 
reemitted at longer wavelengths, whereas most of the longer 
wavelength radiation {A Z 50 fim) just passes through the core 
without interacting at all. The UV, optical and NIR radiation 
that is absorbed is mainly responsible for the heating of the 
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Fig. 7. Density profiles (a), dust temperature profiles (b) and SEDs (c), for BE spheres at T = 10 K under external pressure 
Pext = lO^cm"^ K and with masses 3.5 M0; one subcritical (model BE2, solid lines) and one supercritical (model BE2.2, dashed 
lines). Also for a subcritical BE sphere at T = 15 K with mass 2.6 M0, under external pressure f ext = 10^ cm""* K (model BE7, 
dotted lines). The dash-dot line on the SED graph corresponds to the background SED. The temperature at the centre of more 
centrally condensed cores is lower and the core emission is shifted towards longer wavelengths. 




b (AU) b (MJ) b (AU) 

Fig. 8. Intensity profiles at 90 (a), 170 (b) and 450 yum (c), for the models in Fig.0(BE2: solid lines, BE2.2: dashed lines, BE7: 
dotted lines). The horizontal solid lines on the profiles correspond to the background intensity at each wavelength. At 90 fj.m the 
intensity increases towards the edge of the core but the emission is just 5-10 MJy sr"' above the background and, thus, the cores 
are barely detectable. At 170 fim the intensity drops towards the edge of the core, or rises by a small amount, if the core is cold 
enough. At 450 fim the intensity drops towards the edge of the core in all cases. 



core. A large amount of this radiation will not be available 
if the core is inside a molecular cloud, as we discuss later in 
this paper. The core emits most of its radiation in the FIR and 
submm (also see Fig.Et). The peak of the emission is between 
110 and 160 //m (note that this is the peak of AF^ not F^). 
At these wavelengths, the core is easily observable against the 
background. At shorter wavelengths (e.g. 90 //m) the contrast 
with the background is not very distinct. Finally, in the opti- 
cal the core is seen in absorption and appears like a black blob 
against the bright background. 

The radial intensity profile of a core at a specific wave- 
length /lobs depends on whether this wavelength is shorter or 
longer from the peak wavelength /Ipeak- If ^^obs is much longer 
than /Ipeak (e-g- at 450 iJ.m) then the Rayleigh-Jeans approxima- 
tion for the Planck function holds and the intensity is propor- 
tional to the product of column density with temperature. The 
density decrease towards the edge of the core is much larger 



than the corresponding temperature increase, and the intensity 
decreases considerably towards the edge (see Fig. |8j;). If /Jobs 
is much shorter than /Ipeak (e.g. at 90 fim) then the Wien ap- 
proximation holds and the intensity depends on the tempera- 
ture exponentially, so even a small increase in the temperature 
can balance the density decrease and the intensity increases 
slightly (~ 5 MJysr"') towards the edge of the core (Fig.|S^). 
However, the contrast between core and background radiation 
is very small (~ 5 - 7 MJysr"') and cores should be barely 
detectable at 90 fim. This result is consistent with observations 
of prestellar cores (Ward-Thompson et al. 2002) that show that 
cores are usually well defined at 170 and 200 fj.m but not always 
well defined at 90 //m. Finally if A and /ipeak are comparable 
(e.g. 170 fj.m) the intensity either drops from the centre to the 
edge (A a bit longer than /Ipeak; Fig. 03, models BE2 and BE2.2) 
or it increases (A a bit smaller than /Ipeak; Fig.|8j3, model BE7). 
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Fig. 10. Temperature profiles (a) and intensity profiles at 0.55 fim (b), 170 and 450 //m (c) for a Bonnor-Ebert sphere (^out=l 1-8, 
M - 4 Mq, T - 11 K, fext - 10'*cm"^ K) with different dust scattering properties. The dashed line corresponds to mean 
scattering cosine g = 1 (forward scattering, i.e. in effect, no scattering), the dotted line to g - 0.4 and the solid line to g = . The 
dash-dot horizontal lines on the intensity profiles correspond to the background intensity at the wavelenght noted on the graph. 
Different dust mean scattering cosines do not greatly affect the dust temperature profile in the core. 
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Fig. 11. Temperature profiles (a) and intensity profiles at 0.55 fjm (b), 170 and 450 fim (c) for a Bonnor-Ebert sphere (^out=l 1-8, 
M - 4 Mq, r = 11 K, f ext = 10'*cm"^ K) with different dust scattering opacities. The dashed line corresponds to zero scattering 
opacity, the dotted line to half the MRN scattering opacity and the solid line to the MRN standard model scattering opacity. The 
dash-dot horizontal lines on the intensity profiles correspond to the background intensity at the wavelength noted on the graph. 
Different dust scattering opacities do not greatly affect the dust temperature profile in the core. 



In general, the contrast with the background is quite large at 
these intermediate wavelengths. 

4.3. Effects of Dust Scattering Properties 

The properties of dust in molecular clouds and prestellar cores 
are quite uncertain (see Andre, Ward-Thompson & Barsony 
2000). In this section we examine the effect of different dust 
scattering properties on the temperature profiles and on the 
spectra of prestellar cores. We perform radiative transfer calcu- 
lations using Phaethon, for a supercritical Bonnor-Ebert sphere 
(^out=l 1-8, M = 4 M0, r = 1 1 K, Pext = lO-^cm-^ K) with total 
visual optical depth ry = 30.6, and different dust properties. 

Initially, we vary the mean scattering cosine g. We see 
(Fig. II Oh ) that when the scattering is isotropic {g = 0, solid 
line) the dust temperature at a specific radius inside the core 
is a bit lower (~ 0.3 K) than for the case of forward scatter- 



ing (g - I, dashed line). The case ^ = 1 is equivalent to no 
scattering, so at optical wavelengths there is significant inten- 
sity only at the very edge of the core where the optical depth 
through the core is small and radiation can pass straight through 
(Fig. llOb). If photons are scattered forward, they are able to 
penetrate deeper inside the core and heat it to higher tempera- 
tures. As a result more optical photons are absorbed and more 
FIR photons are emitted. The intensity difference between dust 
models with different mean scattering cosine is very large in the 
optical region (Fig.llOb) but it is only ~ 10 - 20% at FIR and 
submillimetre wavelengths (e.g. 170 and 450 jum, Fig.llOb). 

Next, we vary the scattering opacities of the dust ("Fig. II U . 
We perform 3 calculations: (a) with MRN scattering opacities 
^scat = K^T (solid lines), (b) with k,,,, = (dotted Unes), 

and (c) with no scattering at all (/fscat - 0, dashed lines). The 
results are similar to the previous case: when there is no scat- 
tering (which is the same as g = 1 in Fig.llO> more photons are 
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absorbed by the core, heating it to sHghtly higher temperatures. 
Scattering provides photons with a quick way out of the core 
without them being absorbed. 

This study shows that different dust composition, as re- 
flected in different dust scattering opacity and different scatter- 
ing mean cosine, resuhs in only sHghtly different temperature 
profiles. The optical intensity profiles are strongly dependent 
on the dust scattering properties but at FIR and submillimetre 
wavelengths, where prestellar cores emit most of their radia- 
tion, the intensity is not affected significantly. Thus, we con- 
clude that the scattering properties of the dust do not greatly af- 
fect the results of our radiative transfer calculations of prestel- 
lar cores. 



5. Prestellar Cores Embedded in Molecular Clouds 

In many cases prestellar cores are embedded deep inside 
molecular clouds and the radiation incident on them is differ- 
ent from the interstellar radiation field, and anisotropic (Mathis 
et al. 1983). The ambient molecular cloud acts like a shield 
to UV, visual and NIR interstellar radiation, absorbing and re- 
emitting it in the FIR. It also makes the radiation incident on 
the core anisotropic because in general the molecular cloud is 
not homogeneous and it is not spherically symmetric. Even for 
a spherical ambient molecular cloud with uniform density the 
radiation incident on the prestellar core is not isotropic, even if 
the core lies at the centre of the molecular cloud. That is be- 
cause there will be more radiation incident on a specific point 
on the embedded core from the radial direction (which is closer 
to the boundary of the cloud) than from the tangent or any other 
direction (see Fig.ll2>. 

Another factor contributing to the anisotropy of the radia- 
tion incident on a prestellar core is the presence of stars or other 
luminosity sources in the vicinity of the core. For example, ac- 
cording to the Liseau et al. (1999) model there is a B2V star 
close to p Ophiuchi that increases the UV radiation incident on 
the cloud from one side. Also the NGC 2068/2071 protoclus- 
ters in Orion B (Motte et al. 2001) are in an environment rich in 
FIR, submm and mm photons, from reprocessed UV radiation 
from the newly born stars in Orion. In such cases, the BISRF 
is probably not a very good representation of the radiation field 
incident on the core. 

Previous studies (Evans et al. 2001, Young et al. 2002) have 
acknowledged that deviations from the BISRF are important 
and have used a scaled version of the BISRF that is either en- 
hanced at all wavelengths or selectively at UV and FIR. This 
simple approach has a free parameter, the ISRF scaling fac- 
tor, that is varied arbitrarily to fit the observations but it is not 
connected directly to the molecular cloud in which the core is 
embedded or the transport of radiation inside the cloud, and 
does not account for the fact that the radiation field incident on 
an embedded core is not isotropic. Here, we present more con- 
sistent models in which deviation from the BISRF is a direct 
result of the presence of a molecular cloud that surrounds the 
core. 




Fig. 12. Schematic representation of a prestellar core embedded 
in a molecular cloud (not in scale). The radiation incident on 
the core is not isotropic because rt > Tr = rdoud- 

5.1. Model Description 

We examine the simple model of a spherical prestellar core that 
is at the centre of a spherical molecular cloud (see Fig.ll2>. We 
try to mimic the conditions in p Ophiuchi, where condensations 
have masses in the range 0.05-3 M0 and dimensions 1-6x10^ 
AU (~ 7 - 42 arcsec). The thermal pressure at the edge of the 
cloud is ~ 10^ cm"-' K and the estimated particle density is 
~ 2 X 10"* cm"-' (Liseau et al. 1999). In our study we examine 
cores with dimensions 4 - 8 x lO-' AU and masses 0.4-1.2 Mq. 
We assume that the molecular cloud outside the core has con- 
stant particle density n(H2) - 0.77 x 10"^ cm"^ (coiTespond- 
ing to «tot = 0.96 X 10^ cm"-' for a gas with mean molecular 
weight fi = 2.3 and hydrogen abundance by mass X = 0.7). We 
also assume that the dust in the molecular cloud has the same 
composition as the dust in the core and, therefore, the same 
opacities. As in the study of non-embedded cores we use the 
Ossenkopf and Henning (1994) opacities (see Section l33l . We 
use a BE sphere density profile for the prestellar cores, with 
fixed ambient pressure at ~ 10^ cm""* K. Thus, the free pa- 
rameters in defining the BE profile are the mass of the sphere 
and the gas temperature. Because n{H2) is specified, the visual 
optical depth Tdoud is the only free parameter for the ambient 
cloud. Tcioud also determines the extent of the cloud. Motte et 
al. (1998) calculate Ay ~ 10 mag for p Ophiuchi, but depend- 
ing on the position of the core in the cloud, the extinction could 
be up to ~ 40 mag. We use visual optical depths 5, 10 and 20. 
The detailed parameters of our models are listed in Tabled 

5.2. Temperature Profiles and Mass Estimates 

The dust temperature profile inside the core depends on the 
optical depth of the molecular cloud in which the core is em- 
bedded, and the density profile of the core. (Additionally, the 
dust opacities are important, but we will not study their in- 
fluence here.) The presence of even a moderately thick cloud 
of Tcioud - 5 around the core, shields the core from UV and 
NIR radiation, resulting in a less steep temperature profile in- 
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Table 2. Embedded Prestellar Cores: Model Parameters 



model ID 


M»(Mo) 


T" (K) 


font 


fic^Ccm 


«b"(cm 3) 




^be" (AU) 


T b 
' cloud 


/^cloud" (AU) 


EMI 
EM 1.05 
EMI. 10 
EM 1.20 


0.8 


15 


4.7 


4.5 xlO^ 


6.7 xlO'' 


87 


6.1 XlO'' 



5 

10 

20 


6.1 xlO^ 
1.5 XlO" 
2.3 XlO" 
4.1 XlO" 


EM2 
EM2.05 
EM2.10 
EM2.20 


0.8 


15 


9.4 


2.4 xlO^ 


6.7x10'' 


227 


5.2 xlO^ 



5 

10 

20 


5.2 XlO' 
1.4x10" 

2.3 XlO" 
4.0 XlO" 


EM3 
EM3.05 
EM3.10 
EM3.20 


0.4 


15 


2.4 


1.4x10"* 


6.7 xlO'' 


36 


5.4 XlO' 



5 

10 

20 


5.4x10' 
1.4x10" 
2.3 XlO" 
4.0 XlO" 


EM4 
EM4.05 
EM4.10 
EM4.20 


1.2 


20 


3.6 


1.9 xlO' 


5.0x10'' 


61 


8.1 XlO' 



5 

10 

20 


8.1 XlO' 
2.0x10" 
3.1 XlO" 
5.5 XlO" 


EM5 

EM5.05 

EM5.10 

EM5.20 


1.2 


20 


15.4 


6.1 xlO^ 


5.0x10'' 


430 


6.2 XlO' 



5 

10 

20 


6.2 XlO' 
1.8 XlO" 
3.0x10" 

5.3 XlO" 



" Embedded core properties: M: mass, T: gas temperature, ^out: ^ parameter of the BE sphere > 6.451 for supercritical spheres), 11^: central 
density, n\,: boundary density, Ty: visual optical depth to the centre of the sphere, Rbe' radius of the sphere. 
Ambient cloud properties: T^i^ud'- visual optical the of the cloud (see Fig. ll2t . Rdoud' cloud radius. 



side the core than in the case of a core that is directly exposed 
to the interstellar radiation field. When there is no surround- 
ing cloud the temperature drops from 16 K at the edge of the 
core to around 6-7 K in the centre (AT * 9 - 10 K, depend- 
ing on the core density), whereas with a = 5 cloud the 
temperature drops from around 11 K to 7 K (Ar » 4 K), 
as seen in Figs. 113b and 115b . Particularly, in the case of a 
deeply embedded core (rdoud = 20) the core is almost isother- 
mal (Ar ; 1.5 K) at around 7-8 K, for a not very centrally 
condensed core (Fig. 113b). whereas for a supercritical core 
(Fig. 115b ) Ar a; 3 K. Our studies show that temperatures in- 
side embedded cores are probably lower than 12 K in cores 
surrounded by even a relatively thin cloud (visual Tdoud = 5), 
which seems to be the case for many of the prestellar cores in 
p Ophiuchi. Previous studies (Motte et al. 1998, Johnstone et 
al. 2000) of prestellar cores in p Oph assumed isothermal dust 
at temperatures from 12 to 20 K, when calculating core masses 
from mm observations. At these wavelengths the dust emission 
is optically thin and, consequently, the observed flux is 

Fa = Bi(rdust) n = BATiu.t) ka NiH2)nmH AQ. (9) 
Hence, the inferred column density is 

fiiUH A£2 Ka BAiTdu^iO 

AO. is the solid angle of the telescope beam for a resolved 
source, or the solid angle of the source if unresolved, N(H2) 
is the hydrogen column density, ka is the mm dust opacity, 
and is the Planck function. At mm wavelengths and tem- 
peratures < 20 K the Rayleigh-Jeans approximation holds, so 



fi,i(7^dust) T^dust- Therefore, the estimated column density, and 
consequently the mass, depends on the observed mm flux, and 
the dust opacity and temperature, 

N{H2)^-^. (11) 

Ka i dust 

Thus, the masses of the prestellar condensations calculated by 
Motte et al. and Johnstone et al., using mm continuum obser- 
vations, may be underestimated by up to a factor of 2, which 
will affect their evaluation of the core mass function in the p 
Oph protocluster, and the inferred stability or instability of the 
observed cores. Detailed modelling for each of the prestellar 
cores, taking into account their environment (i.e. surrounding 
cloud and nearby luminosity sources), is needed to calculate 
their masses with more accuracy. Also, as Motte et al. point 
out, the dust opacities and also the dust-to-gas ratio, introduce 
additional uncertainties in mass calculations. 

5.3. SEDs and Intensity Profiles 

In the UV and optical (0.01-1 jum), the radiation coming from 
the system is scattered light and direct background radiation 
(mainly coming from the edge of the cloud where the optical 
depth is smafl). In the NIR and MIR (1-50 yum) most of the ra- 
diation is direct background radiation that just passes through 
the outer, optically thin parts of the cloud. This depends on 
the assumed background radiation field, the optical depth of 
the cloud (and hence the dust properties) and the extent of the 
cloud. The FIR and mm range (60-1300 /im) is the most inter- 
esting area since the core emits most of its radiation at these 
wavelengths. Many terrestrial and space-borne observatories 
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Fig. 13. Density profiles (a), dust temperature profiles (b) and SEDs (c), for BE spheres at T-15 K with mass 0.8 M0, under 
external pressure Pext = lO^cm"^ K, surrounded by a spherical ambient cloud with visual optical depth 20 (model EMI .20, solid 
lines), 5 (model EMI. 05, dotted lines) and (model EMI, dashed lines; no surrounding cloud). The dash-dot line on the SED 
graph corresponds to the background SED. The deeper the core is embedded the lower the dust temperature inside the core. The 
dust temperature is lower than 12 K even when the core is embedded in a relatively thin molecular cloud with visual extinction 5 
mag. 




Fig. 14. Intensity profiles at 90 (a), 170 (b), 450 and 850 yum (c), for the models in Fig.[T3l(EM1.20: solid hnes, EM1.05: dotted 
lines, EMI: dashed lines). The horizontal solid lines on the profiles correspond to the background intensity at the wavelength 
marked on the graph. At 90 /urn the core is seen in absorption against the background but the core is not easily distinguishable. At 
170 fim the intensity increases towards the edge of the core only if the core is not very deeply embedded. However, very sensitive 
observations are needed to detect this feature. At 450 and 850 jum the intensity drops towards the edge of the core. 



cover (or have covered) this range: ISOPHOT//50 (90, 170 and 
200 yum), SCUBA/JMCT (350-1300 fim), IRAM (1300 fim) 
and finally the upcoming SIRTF (3.6-160 /jm, to be launched 
in 2003) and Herschel (75-500 yum, to be launched in 2007). 

At 90 microns the core is seen in absorption against the 
background (Figs. [T4landll6h'). The intensity depends on the 
temperature exponentially, so the relatively small increase in 
temperature towards the edge of the core can compensate for 
the rapid decrease in the density (p ~ r"^), and the intensity 
increases towards the edge of the core. For a very centrally- 
condensed core (e.g. models EM2x, Fig.ll6h) the decrease to- 
wards the centre is around ~ 8 - 10 MJysr ' (depending on 
how deep the core is embedded in the cloud; for deeper embed- 
ded cores the intensity decrease is smaller), and this would be 
very difficult to detect. For less centrally-condensed cores (e.g. 
EMlx, Fig.ll4h') the decrease is even smaller (~ 4-6 MJy sr"'). 



Thus, very sensitive (say ~ 1 - 3 MJysr" ) observations are 
needed to detect cores in absorption. This sensitivity is very 
close to the limits of current instruments, so it is very difficult 
to observe embedded prestellar cores at 90 //m. 

At wavelengths near the peak of the emission (150-250 jum) 
the intensity increases by a small amount (~ 5 - 20 MJy sr"' 
above the background) towards the edge of the cloud and then 
decreases to the background intensity (Figs. [T4b andll6b). If 
the temperature increase towards the edge of the core is big 
enough to compensate for the decrease in density, the outer 
parts of the core are just visible (e.g. models with visual optical 
depth Tcioud = 5). On the other hand, if the increase of the tem- 
perature is not high enough, as happens when the core is deeply 
embedded (rdoud - 20), then the core cannot be distinguished 
from the background. Thus, our models indicate that cores can 
be observed at 150-250 /im only if they are surrounded by a 
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Fig. 15. The same as in Fig.[2J but for a supercritical BE sphere with the same parameters (EM2x models): Density profiles (a), 
dust temperature profiles (b) and SEDs (c), for BE spheres at 15 K, with mass 0.8 M0, under external pressure Pext = lO^cm"^ K, 
suiTounded by a spherical ambient cloud with visual optical depth 20 (model EM2.20, solid lines), 5 (model EM2.05, dotted lines) 
and (model EM2, dashed lines; no suiTounding cloud). The dash-dot line on the SED graph coiTesponds to the background 
SED. 




Fig. 16. Intensity profiles at 90 (a), 170 (b), 450 and 850 yum (c), for the models in Fig.[T5l(EM2.20: solid lines, EM2.05: dotted 
lines, EM2: dashed lines). The horizontal solid lines on the profiles correspond to the background intensity at the wavelength 
marked on the graph. In this case (more centrally condensed core than that in Fig. \l4i . the intensity at the centre of the core 
at 90 i-im is lower and, thus, the core is relatively more easily observed in absorption than a less centrally condensed core. In 
addition, the increase of the intensity towards the edge of the core at 170 fim, is larger in this case and thus easier to observe. 



cloud with a relatively small visual optical depth Tdoud ~ 5, in 
which case the intensity increase is ~ 10 MJy sr"'. Cores could 
in principle be observed even if they are deeply embedded, pro- 
vided there were accurate observations of ~ 1 MJy sr"' at 150- 
250 jum. This result agree with the fact that ISO did not detect 
the prestellar condensations in p Oph (Andre, Ward-Thompson 
& Barsony 2000). 

Finally, at submillimetre and millimetre wavelengths (400- 
1300 fim) the Rayleigh- Jeans approximation for the Planck 
function holds, and the observed intensity is proportional to 
the product of the core column density and temperature. Thus, 
at the edge of the core the intensity drops considerably be- 
cause the temperature increase cannot compensate for the den- 
sity decrease. The core can be easily observed at 400-500 /urn, 
where the contrast with the background is quite considerable 
(~ 50 - 150 MJy sr"'). At wavelengths longer than ~ 600 yum 
the background radiation becomes important and the core emis- 



sion is not much larger than the background emission. For ex- 
ample at 850 jum (Figs.ll4landll6b) the core emission is only ~ 
20-50 MJysr ' above the background, depending on the den- 
sity profile of the core and how deeply the core is embed- 
ded inside the molecular cloud. High accuracy observations 
are needed to observe cores at mm wavelengths, but they are 
available. For example, the sensitivity of IRAM is around ~ 1 
MJysr"'. The peak luminosities at 1300 fj.m that we compute 
with our models are comparable with the observed luminosities 
ofMotte etal. (1998). 

To check if our results are affected by the extent of the am- 
bient cloud, we study a core with the same parameters as the 
EM2.05 model but embedded in a more extended, less dense 
cloud. In both models the optical depth of the cloud is the same. 
As seen in Fig.^J the temperature and intensity profiles at 90, 
170 and 450 jum of the two models are almost identical inside 
the core. This result indicates that the only parameter of the 
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ambient cloud that is important in determining the dust tem- 
perature and the SED of a core embedded in the centre of a 
molecular cloud, is the optical depth of the molecular cloud. 

5.4. Diagnostics 

In Table |3l we list the peak intensities (maximum intensity 
above or below the background) at wavelengths 90, 170, 450, 
850, 1300 yum, for cores embedded in molecular clouds with 
visual optical depths 5 and 20. The lower intensity values cor- 
respond to less condensed cores (subcritical) and the higher 
intensity values to more condensed cores (supercritical). This 
table indicates that embedded cores are most easily distin- 
guished from the background radiation around 450 fim. The 
peak emission from embedded cores could be as low as ~ 10 
MJy sr"' above the background at 1300 fim, but it's at least ~ 
40 MJy sr"' at 450 fim. The wavelength range between 400- 
500 seems favourable for observing embedded cores but the 
atmospheric transmission is not good in this range and space 
observations are needed. The upcoming Herschel space tele- 
scope will be operating in this range. 

Continuum intensity observations at different submm 
wavelengths can be used to determine if a core is subcritical or 
supercritical, assuming cores can be described as BE spheres. 
If the core is supercritical, it is more condensed in the centre 
and the optical depth to the centre of the core is larger than 
a subcritical, less condensed core. This means that the dust 
temperature at the centre of the core is less for a supercriti- 
cal core, and the resultant spectrum is shifted towards longer 
wavelengths. Thus, a supercritical core emits more radiation at 
longer wavelengths (e.g. 450 ^m; see Figll8> than a subcritical 
core. We can exploit the fact that the intensity at 170-200 /^m 
varies little for different cores in the same environment, and 
use the colour index Cl=m(450//m) - TO(170;um) to distin- 
guish between subcritical and supercritical cores, eliminating 
this way any uncertainties about the distance of the observed 
cores. The CI will be larger for supercritical cores. This result 
can be used to determine whether a core is subcritical or super- 
critical when the core cannot be resolved and the usual density 
criterion (pcentre/Pedge > 14.1 for supercritical spheres) is not 
useful. 

Our models can also be used to estimate the visual extinc- 
tion of the ambient cloud surrounding an embedded core. The 
outer parts of even deeply embedded cores (Tdoud ~ 20 - 30) 
are expected to be just visible in emission at 170-200 ^m. (~ 3 
MJysr ' above the background), whereas cores embedded in 
a moderate-thick cloud (rdoud ^ 5 - 7) will be more visible 
(~ 10 MJysr"' above the background), as seen in Figs. 114b 
and 1 16b. The higher the increase in the intensity near the core 
boundary, the less embedded is the core. Thus, very sensi- 
tive observations of embedded prestellar cores at 170-200 //m 
(~ 1 - 3 MJy sr '), might allow us to determine the extinction 
of the cloud surrounding the core, and thus to estimate roughly 
the position of the core inside the molecular cloud. However, 
more sophisticated modelling is required, with a more detailed 
density profile for the cloud and taking into account the close 
environment of the core under study. 



Table 3. Typical peak* intensities for embedded cores 





// (MJy sr-') 




Tcloud = 5 Tcloud = 20 


901° 


5-15 ~3 


170 


10-15 ~ 3 


450 


55-160 40-130 


850 


20-80 15-70 


1300 


10-40 10-25 



* The term peak refers to the maximum intensity above or below 
the background (as noted) at a specific wavelength. 

" These are typical approximate peak intensities for a core embed- 
ded in a cloud with visual optical depth t = 5 and t = 20. The 
deeper the core is embedded the less distinct from the background 
is. The lower value corresponds to a subcritical core and the higher 
value to a supercritical core. 

At 90 yum the core seen in absorption against the background. 



6. Summary 

We have implemented a Monte Carlo radiative transfer method 
to study non-embedded and embedded prestellar cores. This 
method discretises the radiation of one or more sources with 
a large number of monochromatic luminosity packets that are 
injected into the system and interact stochastically with it. Our 
code has been tested against benchmark calculations of other 
well established radiative transfer codes. We have also tested 
our code against the thermodynamic equilibrium test (the sys- 
tem is illuminated by an isotropic blackbody radiation field and 
acquires the same temperature as it) and the pure scattering test 
(only scattering interactions are allowed and the output inten- 
sity is the same as the input one). 

We studied cores that are directly exposed to the ISRF 
and found similar results (temperature and intensity profiles) to 
Evans et al. (2001), using a different radiative transfer method. 
We extended our study to cores that are embedded inside spher- 
ical molecular clouds. We assumed that the ambient cloud has 
uniform density, and that the dust composition is the same as 
that in the embedded core. In this case, the radiation incident on 
the embedded core is not isotropic and cannot be represented 
by the Black (1994) approximation, since the ambient cloud 
shields the core from UV, optical and NIR photons and en- 
hances the FIR and mm part of the spectrum. We found that, 
in this case, the temperature is generally less than 12 K, even 
for an ambient cloud with low visual extinction (~5 mag). The 
temperature gradients inside an embedded core are smaller than 
in the case of a non-embedded core; deeply embedded cores are 
almost isothermal. Recent studies (Andre et al. 2003) using a 
different approach, in which they estimate the efifective radia- 
tion field incident on an embedded core from observations, also 
find that the temperature inside embedded cores is lower than in 
non-embedded cores. Previous mass estimates using mm fluxes 
have assumed isothermal cores at temperatures 12-20 K and, 
consequently, they may have underestimated the masses of the 
cores by up to a factor of 2. However, a more detailed mod- 
elling is needed for each specific core for more accurate mass 
estimate. 
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Fig. 17. Density profiles (a), dust temperature profiles (b) and intensity profiles (c) at 90, 170 and 450 jj.m (the horizontal solid 
lines on the profiles correspond to the background intensity at the 170, 90 and 450 /im, from top to bottom), for a supercritical 
BE sphere at 15 K, with mass 0.8 M0, under external pressure Pext = lO^cm"^ K, surrounded by a spherical ambient cloud with 
visual optical depth 5 and ntot = 0.96 x 10"* cm~^ (model EM2.05, solid lines) and «tot = 0.45 X 10"* cm"^ (same optical depth as 
before but more extended cloud, dotted lines). The dust temperature and the intensity profiles are almost identical inside the core 
for the two models examined, indicating that these profiles are determined mainly by the optical depth of the core rather than its 
physical extent. 




r (AU) r (AU) b (AU) 



Fig. 18. Comparison of a subcritical (model EMI. 05, dotted lines) with a supercritical (model EM2.05, solid lines) BE sphere, 
embedded in a molecular cloud with visual optical depth 5 . Density profiles (a), dust temperature profiles (b) and intensity profiles 
(c) at 90, 170 and 450 /im. The horizontal solid lines on the profiles correspond to the background intensity at the 170, 90 and 450 
fim, from top to bottom. The dust temperature inside a more centrally condensed core is lower than for a less centrally condensed 
core. The core emission is shifted towards longer wavelengths and, thus, a supercritical core will emit more radiation at submm 
wavelengths than a subcritical core does. 



Our models provide a view of cores at a wide range of 
wavelengths. We found that the best wavelength range to ob- 
serve embedded cores is 400-500 yum, where the core is easily 
distinguished from the background. Embedded cores could also 
be observed at 600-1300 fim. The contrast of the core radiation 
against the background radiation is not large but very sensitive 
observations are available in this range. At shorter wavelengths 
the cores are just visible in emission (170-200 //m) or in ab- 
sorption against the background. We also found that very sen- 
sitive observations at 170-200 /im, could be used to estimate 
the visual extinction of the cloud surrounding a core, and thus 
get a rough idea of where the core lies in the in the environ- 
ment of the protocluster. The upcoming Herschel satellite wiU 



be observing in the 60-700 jum range with high sensitivity and 
high angular resolution (Andre 2002), and will test our models. 
Sensitive intensity observations in this range will also reveal 
very low-mass condensations present in embedded protoclus- 
ters, that were previously undetected or poorly detected. These 
observations combined with theoretical models will enable us 
to estimate with great accuracy the temperature profile of re- 
solved prestellar cores. In addition, mm observations from the 
ground will provide accurate mass estimates for the cores, and 
they will constrain the dust opacity. 

Theoretical modelUng should be done for each core indi- 
vidually, taking into account the core surroundings (ambient 
cloud, local luminosity sources). The Monte Carlo approach 
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for the radiative transfer that we have implemented here is in- 
herently 3D and can treat such asymmetric systems. We plan 
to extend our study to ellipsoidal cores, asymmetric ambient 
clouds and anisotropic illuminating radiation fields and, hence, 
model specific cores in embedded protoclusters. 
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